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ABSTRACT 

We study the evolution of the isT-band luminosity function (LF) and the halo occupa- 
tion distribution (HOD) using Subaru observations of 15 X-ray clusters at z = 0.8 — 1.5 
and compare the results with mock clusters (0 < z < 1.3) extracted from the Millen- 
nium Simulation and populated with galaxies by means of the semi-analytic model 
(SAM) of Bower et al., matched in mass to our observed sample. By fixing the faint-end 
slope (a = —0.9), we find that the characteristic luminosity K* defined by a Shechter 
LF is consistent with the predictions of the SAM, which are found, for the first time, 
to mimic well the evolution of A"* in rich clusters at z ^ 1. However, we cannot dis- 
tinguish between this model and a simple stellar population synthesis model invoking 
passive evolution with a formation redshift (zf ~ 5) - consistent with the presence 
of an old red galaxy population ubiquitous in rich clusters at z — 1.5. We also see a 
small difference (AK* ~ 0.5) between our clusters and studies of the field population 
at similar redshifts, which suggests only a weak dependence of the luminous (L L*) 
part of the LF on cluster environment. Turning to our study of the HOD, we find that 
within a radius corresponding to a density 500 times critical, high-z clusters tend to 
host smaller numbers of galaxies to a magnitude K* + 2 compared to their low- z coun- 
terparts. This behavior is also seen in the mock samples and is relatively insensitive 
to the average mass of the cluster haloes. In particular, we find significant correlations 
of the observed number of member cluster galaxies (N) with both z and cluster mass: 
N(M, z) = (53±l)(l + z)" - 61 -" -.2o(m/io 14 -3)0.86±o.05^ p ma u y) we examine the spatial 
distribution of galaxies and provide a new estimate of the concentration parameter 
for clusters at high redshift (c g = 2.8i 'g). Our result is consistent with predictions 
from both our SAM mock clusters and predictions of dark matter haloes from the 
literature. The mock sample predictions rise slowly with decreasing redshift reaching 
Cg = 6.3±°;11 at z = 0. 

Key words: cosmology: large scale structure - galaxies: clusters: general - galaxies: 
evolution - galaxies: haloes - galaxies: luminosity function, mass function. 



1 INTRODUCTION 

The cold dark matter (CDM) paradigm is able to predict 
the formation of structures, which form through grav- 
itational instabilities and cluster hierarchically. A way 
to study the structure formation, is by using the halo 
model formalism, according to which all galaxies in the 
Universe dwell in virialised units of mass called haloes 

* E-mail:dc@astro. livjm.ac.uk 



which obey un i versa l scaling relations. N-body simulations 
|Navarro et all |2004| ) show that dark matter haloes seem 
to have universal velocity and density profiles. The same 
cannot be said for galaxies, which are observed to be 
biased tracers of the mass distribution and whose clustering 
amplitude seems to depend on their properties. Studying 
the bias can provide insight into the physics of galaxy 
formation and the relationship between the galaxy and 
dark matter distributions, with the aim of uncovering the 
processes responsible for the way in which haloes, whose 
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properties are specified by the cosmological mo del, are 
populated with galaxies (|Berlind fc Weinberg||20o3 ). 

A fundamental tool for investigating the influence of 
bias on galaxy clustering statistics is the halo occupation 
distribution (HOD), generally used to determine the average 
number of galaxies within a dark matter halo as a function 
of the halo mass. Within the HOD framework, the virialised 
dark matter haloes with typical overdensities of A ~ 200 
(defined with respect to the critical density) are expected to 
be in approximate dynamical equilibrium and are described 
in one of three ways: (i) P(N|M), the probability of a halo 
of given mass M having an average number of galaxies 
(N); (ii) the relation between galaxy and dark matter 
spatial distributions within haloes and (iii) the relation 
between galaxy and dark matter velo city distributions 
withi n haloes. Theor et ical s tu dies (e.g.. IPeacock fc Sniitr 



when the Universe w as denser (Navarro. Frenk fc White 
19971 hereafter NFW iNeto et al] 120071 ; iGao et all 120081 ; 



Munoz-Cuartas et alj|201ll ). For this reason, understanding 



200d; [Benson et all |2000|; IBerlind fc Weinberg! 



the dependence of concentration on mass and redshift and 
how the relationship between the concentration of dark 
matter and galaxies within haloes evolves with cosmic time, 
can provide fundamental insight into the formation and 
evolution of structures and into galaxy formation process. 

Several statistical galaxy properties have been used to 
empirically measure the HOD assuming a power-law form 
({N) M oc M 13 ) for satellite galaxies. For instance, studies 
of the HOD have been carried out by using either the lu - 
minosity function (e.g., lYang. Mo fc van den Boschl 120081). 
the s patial clustering (e.g., Phleps et al. 20061 ; Abbas et alj 
|2010|) or by counting galaxies within k nown dark mat- 
ter haloes, such as rich clusters (e.g., [k ochanc k et alj 
20031; iLin. Mqhr fc Stanford! 12004 ICollister fc Lahavl 120051 ; 



Berlind et al.ll2003l ; iKravtsov et alj|2004 ) based on N-bodv Ipopesso et al J 120071 ; iHo et ; alj 12009 ). Our study focuses 



simulations, hydrodynamic simulations and semi-analytic 
models (SAMs) showed that P(N|M) can be well modeled 
by a Poissonian distribution in the high halo mass regime 
(Mi, > 10 12 M Q ) while it assumes significant sub-Poissonian 
behavior, which can be modeled by a "nearest integer Q 
(Nint) distribution instead, for lower halo masses. These 
studies also showed (N}m presents a sharp cutoff at low 
halo mass, a slowly rising plateau and a steep increase of 
the occupation number at high halo mass. A simple way 
to model the complicated shape of {N)m, avoiding the use 
of models involving a large number of parameters, is by 
assuming the existence of two separate galaxy populations 
within haloes: (i) central galaxies and (ii) satellite galaxies. 
This choice i s motivated by reas ons based on hydrodynamic 
simulations (|Berlind et al.ll2003l ) and on studies of observed 
galaxy clusters and groups, which take the brightest cluster 
galaxies (BCGs) as a different population from the rest of 
the cluster galax ies. These two popu l ations can be modeled 



separately (e.g., Benson et al. 1 12000| ; IKravtsov etafl [2004: 



I Zheng. Coil fc Zehavi 2007 ). with the simplest case being 
modeling the HOD of central galaxies as a step function 
(N c ) = 1 above a halo mass limit and a power law for 
satellite galaxies. In this way the HOD becomes a measure 
of the combined probability that a halo of mass M hosts a 
central galaxy and a given number iV s of satellite galaxies. 
Studies have also been carried out on the radial, spatial 
and velocity distributions of galaxies and dark matter, the 
majority of them focused on the latter. For instance, the 
halo concentration (defined as one of the shape parameters 
of the radial density profile) and its relation with mass 
and redshift hav e been the topic of many investigations in 
recent years (e.glNagai fc Kr avtsov 20 flj; INeto et al.|[2007l ; 
IGao et alj 12008 ; iMunoz-Cuartas et alj 1201 j ). The impor- 
tance of the concentration resides in its connection with 
the mean density of the Universe at the time of collapse, 
i.e. more concentrated structures formed earlier in the past, 



Referred to as "Aver age" in iBenson et al. and 
IBerlind fc Weinberg! hOOtf ). The definition of this distribution is 
p(Ni\(N}) = 1 - {(N)-Ni), p(Ni + l\(N)) = (N)-N u where N x 
is the integer satisfying N\ < (N) < N\ + 1, with p(N\{N))=0 for 
all other values of N. 



this latter method, which is perhaps the most direct, using 
a sample of 15 X-ray selected clusters at high redshift. 

Unfortunately, there is a significant disagreement 
among the results so far obtained and firm conclu- 
sions are hard to draw. For example, the slope of the 
N — M relation is seen ranging between ft = 0.55 
llMarinoni fc Hudson! 120021 ) and f) ~ 1.7 jAbbas et all 
l20ld) ; while the average concentration parameter of the 
galaxy density profile (c g ) is found to r ange between c E ~ 2 
and c E ~ 8 ([Carlberg: et alj 19971; Ivan der Marel et al 



20001 ; 
2004; 



Biviano fc Girardil 120051; [Katgert. Biviano fc Mazure 



Lin. Mohr fc Stanford! 12004 iRines fc Diaferid [2006 



Muzzin et alj l2007al ; iBiviano fc Poggiantil |2010| ). Although 
the majority of studies have been conducted at low redshift, 
these problems of consistency among different studies of 
the HOD extend to intermediate redshift s (e.g., ILin et all 
120061 ; iBuote et al1l2007l ; lAbbas et aflbOloT ), therefore there 
is a continuing motivation to investigate the evolution of 
the HOD with new galaxy samples. 

The _K"-band LF can be used as a surrogate of the galaxy 
mass function because it is a sensitive probe of the bulk 
properties of galaxy populations out to z ~ 1.5. There 
are several advantages to using the near-infrared light for 
such studies: i) if-band (2.2 /im) luminosities broadly re- 
flect the total stellar mass of galaxies, resulting in a M/L 
ratio that is insensitive to the star-formation history of 
early-type galaxies; ii) stars are easy to remove, as they 
generally have J — K < 1 (Vega), whilst the k-correction 
makes the observ ed colours of the great majority of galax- 
ies J - K > 1 (|De Propris et alj Il999l ; iMcCracken et alj 
|2010| ); iii) k-correction in the near-infrared ba nds varies 
slowly with z and depends wea kly on Hubble type l|Poggiantil 
ll997l ; lBruzual fc Charlotll2003l) ; iv) the effect of extinction at 
these wavelengths is significantly smaller than in optical and 
UV passbands. With the advent of infrared surveys like the 
UKIR T Infrared Deep Sky Survey (UKIDSS. ILawrence etahl 
l2007h and Spitzer, it has been possible to carry out stud- 
ies of the stellar mass of high-2 galaxies. However, such 
studies produced res ults in contrast with th e prediction of 
SAMs . For instance. |Pe Propris et alj (ll999l) ; lEllis fc Jones! 
(|2004l );l Strazzullo et al. I (|2006l) and lLin et al j (feoOrj ) studied 
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the evolution of the cut-off magnitude (K*) of the cluster 
if-band LF out to z ~ 1 and found agreement with passive 
evolving models which have formation redshift (zt ~ 2 — 5), 
suggesting that the bulk stellar mass of K* cluster galax- 
ies has not increased substantially since z = 1. Several 
other studies showed that the high-mass end of the galaxy 
mass function seems to remain pr etty much unchanged since 
z ~ 1 for elliptical galaxies (e .g., ICimatti. Daddi fc Renzinil 
l200rj and |Pozzetti et al.ll20ld) and similar results have been 
obta ined for BCGs (IWhilev et all |200S| ; ICollins et all |2009| 
and IStott et~aL I l2O10f T" These results suggest a timescale 
for the mass assemblage of galaxies similar to the age of 
their component stars, consistent with a monolithic-like 
model, and such activity can be viewed at le ast as qua! 



itatively as consistent w i th a "downsizing" llCowie et al 
' Thomas et al.ll2005l; |Pe Lucia et al.ll2007l ; IStott et a! 



199£ 




2007 





Capozzi. Collins fc StottJl201CI and references therein) 



process, according to which the more massive early-type 
galaxies end their star formation and settle on the colour- 
magnitude relation earlier than their less massive counter- 
parts. 

By contrast, the SAM used by |Pe Lucia et all (120061 ) 
predicts that the majority (~ 70 per cent) of the stellar 
mass of ellipticals at z = is already formed by z = 1, 
but that at this redshift only a few per cent of this mass is 
assembled in the main progenitor. 

In this paper we investigate the process of galaxy 
formation by studying the if-band LF (evolution of K*) 
and the HOD of a sample of 15 galaxy clusters, containing 
the majority of the highest-z (z > 0.8) X-ray clusters 
observed so far. We want to push the study of the evolution 
of cluster galaxies to higher z by extending the Hubble 
diagram of K* out to z ~ 1.5, because it is at z > 0.8 that 
the differences among the evolutionary predictions based on 
stellar populations models are enhanced. The HOD of this 
sample is investigated using the number (JVsoo) of cluster 
galaxies within -R500 (the radial distance where A = 500) 
as a function of cluster mass and z. We also study galaxy 
clustering by analysing the radial galaxy surface number 
density profile in order to estimate the galaxy concentration 
parameter at z ~ 1 to compare with sim ilar estimates 
at low z (e.g., iLin. Mohr fc Stanford! l2004h ■ We compare 
the results of our HOD analysis with results obtained 
from dark matter haloes of similar redshift selected from 
the Millennium Simulation (MS; ISpringel et al. 

I l2005h . 

whose hal oes are popu l ated w ith galaxies taken from the 
SAM by iBower etall (|2006h . It is worth pointing out 
that the majority of previo us observation-based studies 
estimated c g at low z (e.g, ILin. Mohr fc Stanford! |2004| ; 
iPopesso et al. 2007h. apart from some recent exceptions 
(IBiviano fc Poggiantil |201Ch . In addition, even though 
more attention to the concentration parameter has been 
given by theoretical s t udies based on N-body simulations 
(e.g. iNeto et all [20071; iGao et"afl 120081 ; iDuffv et all 120081 ) 
and SAM (|Nagai fc Kravtsovll2005l ~ the majority of them 
focused more on Cd m (the concentration parameter of dark 
matter haloes' radial density profiles) and its relation with 
halo mass and z. Here, instead, simulations and their 
associated SAMs are used to study the concentration of 
galaxies c g . To our knowledge this is the first attempt to 
carry out a self-consistent comparison between the observed 



galaxy concentration parameter in clusters and that in 
mock clusters within the MS. 

The paper is organised as follows: in Section [5] we 
describe the observed cluster sample, while Section [3] is 
dedicated to data reduction and photometry. Sections [3] 
and [S] arc focused on the study of the if-band LF and the 
HOD respectively of the observed and mock samples, while 
in Sections [5] and we present and discuss our results. 
Finally, we draw our conclusions in Section [H] 

Throughout this paper we make use of magnitudes in 
the Vega photometric system and assume a standard cosmol- 
ogy with H = 70 km s" 1 Mpc" 1 , fi m = 0.3 and fi A = 0.7. 



2 CLUSTER SAMPLE 

The clus ters are a subsam ple of the cluster sample de- 
scribed in IStott et al.1 (|2010h and consist of the clusters from 
that list with av ailable J and K band data. Two clusters 
(CL J1226+3332 IMaughan et al.ll2004 and MS1054.4-0321, 
iBranchesi et al.l 12007 ) were included in the original sam- 
ple, but we decided to exclude them because the MOIRCS 
field of view (FOV) (4' x 7') (see below) was not extended 
enough to contain their large -Rsoo- Our final sample con- 
sists of the 15 clusters between 0.8 < z < 1.5 detailed in 
Table [T] Some of these were discovered by various X-ray 
surveys whilst some were optically selected clusters showing 
extended X-ray emission. All 15 clusters have spectroscopi- 
cally confirmed redshifts (Fig. [1} and X-ray luminosities in 
the range 1 < £x < 19 x 10 44 erg s _1 . Cluste r mass esti- 
mates are made using the M — Tx relation in IStott et~al 



2010), whose parameter values are based on the IMaughan 



1 2007 ) derived M — Tx relation. We refer to the study by 



IStott et al.l (|2010l ) for further details about the derivation of 
cluster masses. 



3 OBSERVATIONS, DATA REDUCTION AND 
PHOTOMETRY 

The observations we re taken with the MOIRCS camera 
(|lchikawa et aljfeood ) on the 8.2 m Subaru telescope, which 
provides imaging and low-resolution spectroscopy over a to- 
tal FOV of 4' x 7' with a pixels scale of 0.117" per pixel. Ob- 
servations were taken in 0.5" seeing on the nights of August 
8th and 9th 2007 and in 0.3" - 0.6" seeing on the nights of 
December 16th 2008 and April 18th 2009, with the clusters 
centered on Detector 2. A circular 11-point dither pattern 
of radius 25" was used for both bands to ensure good sky- 
subtraction. The modal integration times were 25 min at J 
and 21 min at K s , although we observed some of the higher 
redshift clusters for 50 per cent longer when scheduling al- 
lowed. These exposures reach a 5a limiting magnitude of at 
least J=21 and K s = 20 (Vega). 

The data are reduced using the external IRAF package 
MCSRED. They are flat-fielded, sky subtracted, corrected 
for distortion caused by the camera optical design and 
registered to a common pixel coordinate system. The final 
reduced images on which we perform the photometry 
are made by taking the 3cr (s.d.) clipped mean of the 
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to ensure accurate colour determination. This photometry 
is carried out by using SExtractor MAG_APER magnitude 
using a circular aperture of diameter 1". 
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Figure 1. Spectroscopic redshift distribution of the 15 observed 
clusters in our sample. 



The photometry is calibrated to the Vega system using 
a combination of standard star observations and the Two 
Micron All Sky Survey (2MASS, ISkrutskie et~alll2006h and 
UKIDSS catalogues. The typical photometric errors are 0.01 
and 0.08 for the standard star and survey calibrated data, 
respectively. 

In Figure[2]the image of CL J0152. 7-1357 obtained after 
all the described process is shown as an example. 



4 ANALYSIS OF THE OBSERVED CLUSTER 
SAMPLE 

In order to carry out our study of the LF and the HOD, 
both foreground and background objects h ave to be re- 
moved or accounted for . It is well known fe.g-. lLegget j|l992l ; 
|Pe Propris et al.l Il999h that stars generally have observed 
J — K < 1 (Vega), while, because of k-correction, the great 
majority of galaxies, excluding the most local ones, lie in the 
region J — K > 1. Carrying out such a colour cut effectively 
removes contamination by stars. However, there is the pos- 
sibility that faint blue galaxies have observed J — K < 1, so 
for this reason we only perform this J — K < 1 at magni- 
tudes brighter than K = 18, i.e. at least 2 mag brighter than 
the 5<J limit of each cluster. Fainter than this, field galaxy 
counts outnumber those of the stars by more than a factor 
of 10 (|De Propris et al J 1 19991 : iMcCracken et aij|2010h , thus 
at these faint magnitudes the correction for stellar contam- 
ination becomes relatively unimportant. From simulations 
we estimate a contamination of only 3 — 4 stars per cluster 
at K < 18, which remains stable down to a 4cr detection 
limit. 

The correction for foreground and background galax- 
ies is performed statistically using the UKIDSS Ultra Deep 
Survey^} (UDS) Data Release 5 field. This field consists of a 
region of A = 0.77 sq. deg in the northern hemisphere with 
a 5<t depth of 21.5 and 22.5 in K and J bands, respectively. 
For reasons of homogeneity, the same method is used for 
carrying out the star-galaxy separation in the background 
field. 

For all galaxies in this field, sky coordinates, total mag- 
nitudes and aperture magnitudes (used to calculate galaxy 
colours) are extracted from the UDS database. 




2'x2' 




Figure 2. 2' X 2' if-band image of CL J0152. 7-1357 



dither frames. The galaxy photometry is extracted using 
SExtractor (version 2.5) MAG_AUTO magnitude, which 
is found to be within ~ 0.1 mag of the total for extended 
sources (|Martinil 1200 ll ). To calculate the colours of the 
galaxies, we run SExtractor in dual image mode so that 
the K s (hereafter K, for simplicity) band detections extract 
the J band catalogue with identical positions and apertures 



4.1 if-Band LF 

The cluster sample is divided into three redshift bins (me- 
dian 2=0.85, 0.97, 1.23) roughly containing the same num- 
ber of clusters. For each cluster, after the stars are removed, 
we determine the A'-magnitude distribution binned in 0.5- 
mag intervals. The number o f galaxies in each b in JV(, in is 
calculated using the formula (jEllis fc Jonedl2004h : 

Nun = N cl - N back ( ) , (1) 

\ Aback J 

2 A detailed description of the survey can be found at: 
http : //www .nottingham. ac .uk/astronomy/UDS/ 
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Name RA Dec z T x (keV) M 20 o (10 14 M Q ) Reference 



RDCS J1317 + 2911 


13: 


17: 


21 


.70 


+29:11:18.0 


0.81 


— U.o 


2. 


7+2.9 
— 1.3 


Branchesi et al. (2007) 


CL J0152.7 - 1357 


01 


52 


11 


.00 


-13:57:45.0 


0.83 


5.4+g-g 
— u.y 


1 


5 +2.7 


Vikhlinin et al. (20091 


CL J1559.1 + 6353 


15: 


59 


06 


.00 


+63:53:00.0 


0.85 


— l.U 


2. 


's+f| 

— 1.0 


Mauehan et al. (2006) 


CL J1008.7 + 5342 


10 


08 


12 


.00 


+53:42:00.0 


0.87 


3 g+0.8 


2. 


9+1.6 
-1.0 


Maughan et al. (2006) 


CL 1604 + 4304 


16: 


04 


25 


.20 


+43:04:53.0 


0.9 


9 c+l.l 
z -°-0.7 


1. 


9+1.6 
z -0.6 


Lubin. Mulchaev & Postman (2004) 


CL J1429.0 + 4241 


li: 


29 


06 


.40 


+42:41:10.0 


0.92 


6 2+ 15 


5 


5+5.3 
°-2.0 


Mauehan et al. (2006) 


RCS J0439 - 2904 


04: 


39 


38 


.00 


-29:04:55.0 


0.95 


i 5 +0.3 
1 '°-0.2 





5+0.4 
°-0.2 


Hicks et al. (2008) 


2XMM J083026 + 524133 


08: 


30 


25 


.90 


+52:41:33.0 


0.99 


o 9+O.9 


8, 


5+4.1 
°-3.4 


Lamer et al. (2008) 


WARPS J1415.1 + 3612 


li: 


15 


11 


10 


+36:12:03.0 


1.03 


6 2 +0 - 8 
°' z -0.7 


5. 


9+2.9 
-1.9 


Branchesi et al. (2007) 


RDCS J0910 + 5422 


09: 


10 


li 


.90 


+54:22:09.0 


1.11 


6.4H : i 


5. 


q+4.1 
°-2.6 


Balestra et al. (2007) 


RX J1053.7 + 5735 (West) 


10 


53 


39 


.80 


+57:35:18.0 


1.14 




2. 


7+1.4 
'-1.0 


Hashimoto et al. (2004) 


XLSS J022303.0 - 043622 


02. 


23 


03 


.00 


-04:36:22.0 


1.22 


q 5+0.4 
°"°-0.4 


1, 


c+0.9 
S -0.7 


Stott et al. (2010), Bremer et al. (2006) 


RDCS J1252.9 - 2927 


12 


52 


54 


.40 


-29:27:17.0 


1.24 


7 9+0.4 


6 


1+2.3 
-2.4 


Balestra et al. (2007) 


XMMU J2235.3 - 2557 


22 


35 


20 


.60 


-25:57:42.0 


1.39 


8.6±i;| 


7 


7+4.4 
'-3.1 


Rosati et al. (2009) 


XMMXCS J2215.9 - 1738 


22 


15 


58 


.50 


-17:38:03.0 


1.46 


4.i±S:S 


2 


i+1.9 

-■--o.s 


Hilton et al. (2010) 



Table 1. The observed cluster sample. 



where iV c i and A c \ represent the number of cluster galaxies 
and the cluster area respectively; while Aback and ^4back are 
the corresponding numbers for the background field. 

The corresponding error in each bin interval is found by 
summing in quadrature the Poissonian error on 7V c i and the 
error on the background - which is made up of a Poissonian 
term and a term accounting for galaxy clustering and given 
by 

°N back = Vn^J 1 + — . (2) 

Here 9 C is the angular radius such that Q, = ty9 c 2 , Q, being 
the solid angle of the background field, and N is the number 
density of galaxies in the bin considered. The error on Aback 
is then normalised to the projected surface area of the clus- 
ter. The parameters S and Au, describe the galaxy angular 
correlation function such that: 

u(6) = A„0 S . (3) 

We adopt t he values 6 = -0 . 8 and A u = (13.49+1.57) x 1CT 4 
taken from lTemporin et al.l (|2008l ). who studied the A-band 
angular correlation function down to K — 20.5, using the 
VIMOS VL T Deep SurveE (VVDS ). 

Similarly to lDe Propris et al.l (1 19991 ). we then use the differ- 
ence in redshift between the cluster redshift and the median 
redshift of its assigned bin, along with the k-correction, to 
appropriately transform the galaxy magnitudes. In this way, 
the intervals will align at the median redshift of each bin. 
The k-correction i s cal culated by using the model of 
iBruzual fc Charlotl (|2003l ). 

For each redshift bin, we obtain a stacked luminosity 
function (Fig. [3} . BCGs are excluded, as they affect the LFs 
in such a way that the abundance of very bright galaxies 
is underesti mated when a Schechter function is used for th e 
fitting (e.g.. ISchech"teilll976l ; Ichristlein fc Zabludofl1l2003h . 
We then fit the stacked LFs obtained in each bin with a 

3 Detailed information about this survey can be found at: 
http : //www . oamp. f r/virmos/vvds . htm 



Schechter function (|Schechterl [l976p . fixing the faint end 
slope to a = —0.9, under the supposition this parameter 
does not evolve with cosmic time and solving for both A* 
and the normaliz a tion. This is the value of a measured by 
|Pe Propris et al.l (|l998h for the A-band LF of the Coma 
Cluster at z ~ 0.02. The choice of fixing a is justified by 
our photometry only reaching ~ 1.5 mag below L*. The fit- 
ting procedure was based on least-squares based on the the 
Levenburg-Marquardt algorithm. The values of the best fit 
A* for the stacked LFs are reported in Table while the 
best fit LFs (plotted as average galaxy number counts per 
cluster vs. magnitude) are shown in Figure [3] 

Since the value o f a measured in the field is a ~ — 1 
(|Cirasuolo et al.ll2010l ) we also repeat the fit for the stacked 
LF fixing the slope to a = — 1. The best fit values of A* 
obtained (Table [2} are consistent within la of the ones cor- 
responding to a = —0.9. The latter values are plotted in the 
Hubble diagram in Figure |4j together with the predictions 
of evoluti onary models with differen t formation z calculated 
using the lBruzual fc Charlotl (|2003l ) model. The models are 
calculated assuming a simple stellar population (SSP) with 
a Chabrier initial mass function and solar metallicity. In 
Figure U we also compare our re sults with the predictions of 
the SAM byH 

owcr et al.l (|2006j ). which implements baryon 
physics, such as AGN feedback from a central black hole 
and star formation, as a SAM bolted onto the MS. Mock 
clusters were selected from the simulations in seven redshift 
snapshots as described in Section[5] The A* values from the 
mock clusters are obtained by applying identical selection 
criteria to the real data, including the same least-squares 
fitting procedure. 



4.2 HOD Analysis 

Despite the diverse methods for studying the HOD, directly 
counting the number of galaxies within known dark matter 
haloes remains the most straightforward way to estimate the 
average galaxy number within a halo as a function of halo 
mass. Furthermore, for galaxy clusters their total mass can 
be estimated via gravitational lensing or, as in this case, us- 
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Figure 3. if-band stacked LFs for our three z bins (z =0.85, 0.97, 1.23) with their best-fitting values of K* . 



ing the temperature of the X-ray emitting gas sitting in the 
cluster's gravitational potential well. Estimating the halo 
masses in s uch a way can constrain qu ite directly P(N|M) at 
high mass l|Berlind fc Weinberdl2002h . Furthermore, the re- 
lation between dark matter and galaxy spatial distributions 
can be investigated in clusters by means of the radial distri- 
bution of their galaxies. The latter can be directly probed 
by their surface density profile. 

As mentioned before in Section [Q several studies (e.g. 
iBerlind fc Weinberdl20o3 ; lKravtsov et alj|2004 ) have shown 
that {N) M can be described as the joint probability that a 
halo of mass M hosts a central galaxy and that the halo hosts 
a given number N s of satellite galaxies. Since this study in- 



vestigates the high mass regime of {N) M , where P(N|(A r }) 
is expected to be Poisson distributed, only satellite cluster 
galaxies are considered; moreover the observed and mock 
clusters analysed here always have M200 > 10 14 Mq where, 
as described in the Introduction, (N) M oc M 13 . In addition 
these clusters (mock and real) are also used to study c g at 
high redshift. We note here that in order to directly com- 
pare our results with those of other workers we transform 
from M200 and -R200 to the respective values at over den- 
sity 500. In this paper the mass transformation is always 
ca rried out using t he st andard NFW profile prescription 
bv lHu fc Kravtsovl (|2003l ) and for the radial scaling we use 
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z 


K* 


K* 


{M 500 > 


(N 5 00> 


No. Clusters 




(a = -0.9) 


(a = -1) 


(1O 14 M ) 


(gal) 




0.85 


17.8±0.3 


17.4±0.6 


2.0 ± 0.4 


48±8 


5 


0.97 


18.3±0.2 


18.0±0.6 


3.6 ± 1.2 


56±7 


4 


1.23 


18.4±0.2 


18.6±0.3 


3.1 ± 0.7 


46±6 


6 



Table 2. Here are tabulated, for each redshift bin, the values of 
median z (col 1), best-fitting K* for a = —0.9 (col 2), best-fitting 
K* for a = —I (col 3), average galaxy number within -R500 (col 
4) and number of clusters per 2 bin (col 5), measured for our 
observed cluster sample. 



A/a = (4/3)7r R^Ap c (z), where A is the considered overden- 
sity. 



4.2.1 N500 - M500 Relation 

After determining the values of K* for each of the three 
redshift bins, we are able to calculate (-/V500) to a depth 
of K — K* + 1.5. However, in order to allow us to mean- 
ingfully compare our results with other studies at low z 
, where at least a depth of K = K* + 2 is reached, we 
decided to extrapolate {N500) to K* + 2 (after normalis- 
ing brighter than K — K* + 1.5 and allowing for suitable 
propagation of errors). The values obtained are shown in 
Table [21 We plot thes e valu es against M500 together with 
iLin. Mohr fe Stanford! (|2004r ) data points and best-fit rela- 
tion at 2 ~ 0.06 in Fig. \E\ 



4-2.2 Galaxy Surface Density Profile 

After selecting only cluster galaxies (BCGs excluded) within 
R/R500 < 1 an d out to K = K* + 1.5, we stack them all to 
obtain a single radial profile. We then calculate the surface 
number density normalised per virial area and plot it as a 
function of virial fraction (i?/i?20o)- In this way we are able 
to obtain a stacked radial number density profile and to 
study the mean concentration of cluster galaxies at z ~ 1. 
In order to estimate c g we fit a NFW density profile to the 
data. The three-dimensional NFW profile has the following 
form: 

P^l = ^ (4) 

Pc (r/r s )(l + r/r a ) 2 ' 

where p c = 3E~o 2 /8nG is the critical density for closure, S c 
is a characteristic density contrast and r s is a scale radius. 
Since we consider a density contrast of A = 200, S c and r s 
take the following form, 



S c = 



(200/3)c 



[ln(l + c) - c/(l + c)] ' 
r s = r 2 oo/c, 



(5) 
(6) 



where c is the concentration parameter of the profile. The 
surface density is then obtained by projecting the three- 
dimensional profile along the li ne of sight from ne gative to 
positive infinity and, following iBartelmannl (|l996t ). can be 
written as 



£(x) = 
where p.. 



2p s 



- 1 

Sc Pc 



/(*) 



(7) 
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• This work 

• Como De Propels et ol. 1998 
o De Propris et ol. 1999 

x Muzzin et ol. 2007 
A Ellis & Jones 2004 individuols 
X Ellis & Jones 2004 composite 

□ Strozzullo et ol. 2006 

□ CIrosuolo et ol. 2007 (Field) 




No evolution 
z,=5 
z,= 2 
z,= 1.5 

SAM cluster prediction 
Field best-fit (Cirosuolo et ol. 2010) 



1.5 



Figure 4. Hubble diagram for K* . The dat a points for the 
cluste r envir o nment are taken from: t his work. |Pe Propris et alj 
(Il998l. I1999T), iMuzzin et all j2007al1 . lEllis fc Jones! d2004h an d 
IStrazzullo et all bOOrjl . From the field: ICirasuolo et alj J2007tl . 
Several evolutionary models are also plotted: no evolution in ab- 
solute if-band magnitude (full line), passive evolution with Zf= 
1.5, 2 and 5 (das hed-triple-dot t ed, d ashed-dotted and dashed line 
respectively) and Bower et al. ( 2006) SAM predictions (green full 
line) obtained by using mock clusters (see text). The grey-shaded 
region shows the typical variation in the evolution of the SAM 
K* when a constant halo cluster mass is adopted. The scatter is 
dominated by Poisson statistics. T he parametrizat i on of the ob- 
served evolution of the field K* bv lCirasuolo et all j2010l) is also 
shown (blue full line). 




(8) 



cr/r2oo and fix) is given by: 



The limited FOV of MOIRGS does not allow us to observe 
galaxies far enough from the cluster to be able to deter- 
mine the background level by investigating regions of the 
sky dominated by field galaxies. For this reason, and be- 
cause we do not have spectroscopic redshifts for our clus- 
ter galaxies, we can not directly background subtract our 
stacked numbe r density profile. Hence we start with the 
methodology of ILin. Mohr fc Stanford! (|2004h and carry out 
a maximum likelihood fit to the projected radial density pro- 
file with the supposition of having a constant background. 
The observed profile is binned very finely in virial fraction 
(R/R200) so that in every bin, the density within the differ- 
ential area dS could be considered constant. By doing this 
we can estimate the probability of obtaining the measured 
galaxy counts within dS for a specific prediction (iV mc ,dei) 
given by the model. This situation is well represented by 
a Poisson probability distribution with mean A = -/V mo d e i- 
Our model is simply given by the sum of the cluster contri- 
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• Millennium Simulation 
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• Lin, Mohr & Stanford 2004 
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1.0 
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Best-fitting Profile 



« 10" 
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c = 2.8 

N =553 

N B « k =527 




10' 

0.001 
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Figure 5. -/V500 against Msoo- Our mean values are reported 
as blue full dots (observed sample), red full dots (mock clus- 
ters with 2 > 0) and green full dot (mock clusters at z = 0). 
Note that the observed data point at z = 1.23 has been slightly 
shifted towards lower masses to make it visible. The data points of 
iLin. Mohr &: Stanfor -plotted as grey full dots— together 

with their best-fitting relati on are shown. We note that, since 
ILin, Mohr &: Stanfor d| j2004h values are measured over clusters 
at z ~ 0.06, only the average value for -/V500 at z = (green full 
dot) can be directly compared with them. Note our simulation- 
based measures are average values measured over the 100 most 
massive clusters in each snapshot (see text). 



bution to the projected density and the contribution of the 
background: 



E — E c ; + Efc ac 



(9) 



where E c i is given by eq.[7]and Eback is a constant. Specifi 
cally, the quantity 

N bin 

- log C = - lo S [ 



■e - A 'Ai"°' 



(10) 



is minimised, where -/Vbm is the number of bins used, N h Si 
is the observed number counts in the differential area dSi 
and Ai = E(xi)d5*i is the number of galaxies predicted by 
the model (cluster galaxies plus background, iV mo dei) in the 
same differential area as a function of p s , c = c g and Eb ac k- 
In the fitting routine we allow the parameters p s , c g and 
Eback to vary. The advantage of this method is that the 
likelihood involves a proper sum of Poisson probabilities 
(in this respect our metho d improves on the one used by 
ILin, Mohr fc Stanf ord 2004) and that, unlike methods based 
on the x 2 ) the results are relatively independent of the bin- 
ning. Figure [6] shows the stacked radial density profile out 
to ~ 0.65i?200 corresponding to the -R500 radius considered 
for each cluster. The determined best fitting galaxy con- 
centration parameter is c g = 2. 8*0 s> with the uncertain- 
ties derived from the relation A\ 2 = — 2Alog£. The value 



Figure 6. Stacked radial surface number density profile for our 
observed cluster sample. Galaxies within -R500 and out to K = 
K* + 1.5 only are used. The best-fitting profile is shown as a 
dashed line. 



obtain e d is consistent with t hose fo und by ICarl bcrg ct al 



19971) . Ivan der Marel et al.1 teOOOh. ILin. Mohr fc Stanford 



2004) and Biviano fc Poggiantil ( 2O10l ) for clusters at low 
and intermediate redshifts (c g = 3.7, 4.2, 2.9 and 4.0 at 
2 ~ 0.3, 0.3, 0.1 and 0.55, respectively). We must point 
out though, that the magnitude limit used in these studies 
does not always coincide with ours. 

The background density obtained amounts to around 
Eback ~ 8 per cent of the peak density of the stacked pro- 
file or ~ 49 per cent of total galaxies (527 out of 1080). We 
tested the reliability of these estimates by analysing the ra- 
dial number density profile out beyond i?2oo- This is done 
by using those clusters in the sample (3 out of 15) whose 
-R200 is small enough to allow us to estimate the background 
out to at least 1.4i?200- For these clusters the background 
level is evaluated by averaging the values in the outermost 
bins (R/R200 > 1) °f the radial profile. The values obtained 
vary between 4 and 12 per cent of the density peak value, 
corresponding to 30 — 50 per cent of the total galaxy number 
out to the same radial distance as used in the stacked profile 
(R ~ 0.65i?20o)- This result is in relatively good agreement 
with the background level estimated with the likelihood. 



5 ANALYSIS OF MOCK CLUSTER SAMPLES 
5.1 Cluster Halo Selection 

The MS3 is a large N-body simulation of the standard 
ACDM cosmogony. It follows 2160 3 particles, each of mass 

4 Detailed information about the simulation can be found 
at: http://www.mpa-garching.mpg.de/millennium/ A dc- 
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Redshift /3 (M500) <N 5 oo) c g c dm Alogc 

(1O 14 M ) (gal) 






0.82 ± 0.01 (3769) 


7.6 ±0.3 


138±1 


6 30+ 39 


3.74 


r, r.q + 0.03 

u,io -0.02 


0.21 


0.76 ± 0.01 (2752) 


5.6 ±0.2 


100±1 


c Qi+0.39 
°- aL -0.35 






0.41 


0.82 ± 0.02 (1953) 


4.3 ±0.1 


84±1 


5 05+ - 41 
°' UO -0.38 






0.62 


0.86 ± 0.02 (1289) 


3.2 ±0.1 


45±1 


4 90+ - 52 


3.48 


1 5 +0 - 0S 


0.83 


0.74 ±0.03 (831) 


2.4 ±0.1 


45±1 


oc+0.45 
J-oO_0.40 






0.99 


0.78 ±0.04 (541) 


2.0 ±0.1 


28±1 


o 7n +0.50 
°- ,u -0.44 


3.66 


00+ 06 
u ' uu -0.05 


1.27 


0.81 ±0.08 (233) 


1.40 ±0.05 


21±1 


3.77t°;{g 







Table 3. Properties of the mock cluster sample. Coll: redshift of the snapshot bin; col2: the slope f) of the -/V500 — M500 relation (HOD 
power-law index) for all the mock clusters in the snapshot with the number of clusters shown in parenthesis - note that all other estimates 
from the mock sample in this table use the 100 most massive cluster haloes at each redshift snapshot; col3: the average mass within the 
radius corresponding to A = 500; col4: the average galaxy number within K500 brighter than K* + 2 (these values do not include the 
systematic correction of 10 per cent discussed in Section 15.26 ; col5: the average galaxy concentrati on paramet e r for the mock clusters; 
col6: the average dark matter concentration parameter as measured from the C(j m (M, 2) relation of lGao et al.l ||2008T ). not available for 
all z; col9: Alogc = log(c g /c<j m ), see text. 



8.6 x 10 s h^ 1 M Q , in a cubic box of 500 h" 1 Mpc on a 
side from redshift 2 — 127 to the present. The cosmological 
parameters were chosen to be consistent with a combined 
analysis of the 2dFGRS (|Colless et al.ll200ll ; iPercival etajj 
l200ll ) and first - year Wilkinson Mic rowave Anisotropy Probe 
(WMAP) data (|Spergel et al.ll2005t l within the concordance 
cosmology. 

By scanning the entire volume of the simulations we ex- 
tract all haloes above a minimum mass of M200 > 10 14 M , 
selected in seven snapshots corresponding to redshifts 2 = 0, 
0.21, 0.41, 0.62, 0.83, 0.99 and 1.27, respectively. 

The purpose of the simulations is to follow the evolution 
of the observed clusters, which do not constitute a statisti- 
cal sample but are highly biased to the most X-ray luminous 
systems and therefore the most massive clusters at z > 0.8. 
The finite volume of the MS means that the number of mas- 
sive clusters is rare, so for the study of the evolution oi K* , 
{N500) and c g , we restrict our analyses to the most massive 
cluster haloes in each snapshot of the simulation, since they 
are thought of as equivalent to t he brightest X-ray s y stems . 
This approach follows that of |Pe Lucia fe Blaizotl (|2007l ) 
who used the 125 most massive clusters in similar snapshot 
intervals from the MS to investigate the evolution of BCGs 
in X-ray selected luminous clusters (see Sect. 6 of their pa- 
per). Selecting the most massive 100 clusters per snapshot 
we can reproduce both the observed redshift and mass dis- 
tributions of our high-redshift clusters reasonably well; for 
the mocks: 0.83 < 2 < 1.27, 1 < M 500 < 7.7 x 10 14 M Q and 
median M500 = 1-8 x 10 14 Mq, and for the observations: 
0.81 < 2 < 1.46, 0.4 < Af 5 oo < 6.2 x 10 14 M Q (next least 
massive cluster is at 0.9 x 10 14 Mq, so only one observed 
cluster lies significantly outside the mock mass range), and 
median M500 = 2 x 10 14 Mq. Details of the mock clusters 
used are given in Table [3] 

One consequence of selecting the most massive haloes 
per snapshot is that the m ass range evolves w i th co smic 
time (as in the analysis of |Pe Lucia fc Blaizotl [20071 ). so 
to investigate whether our results are biased by this mass 



tailed description of the data base can be found in 
iLemson &: Virgo Cons ortium ( 2 006h . 



evolution we also present results using other mock cluster 
samples selected with a constant average mass and a con- 
stant mass range over cosmic time. Results of these analyses 
are reported in Section 16.2.31 The galaxies populating the 
haloes used are extracte d from the SAM ( see also Section 
I4.1|l galaxy catalogue by I Bower et al.l (|2006l ). 



5.2 Satellite Galaxies' Power-law Index & 
Occupation Number 

Once the galaxy spatial coordinates and if-band magnitude 
are extracted for the totality of haloes contained in each red- 
shift bin, we select galaxies with projected radial distances 
within -R500 for each mock catalogue and carry out individ- 
ual fits to obtain individual LFs in the same way as described 
for the observations in Section l4.1l We point out that in or- 
der to reproduce observations the cut in -R500 is not applied 
along the line of sight and that galaxies are assigned to each 
halo by using the unique halo IDs, which allow us to univo- 
cally associate galaxies to the haloes they belong to. How- 
ever, on larger scales the lack of a background in the mocks 
means that the galaxy selection differs from the observed 
clusters, for which galaxies are counted within a cylinder of 
projected radius -R500 and a background correction applied 
using galaxy counts (see Section 14. 1[) . In principle line-of 
sight contamination from galaxy clustering on scales extend- 
ing beyond the virial radius of the clusters could affect our 
richness and concentration values by introducing an offset 
between measured and mock estimates. Therefore to check 
this we inspect the number of non-cluster galaxies at 2 = 1 
brighter than K* + 2, with physical and projected distances 
from mock halo central galaxies of -R500 < R < 3 Mpc and 
R < Rboo respectively. We find that -/V500 values for the 
mock clusters are ~ 10 per cent lower on average compared 
to those observed due to this contamination. This is less 
than the individual Poisson error on N500 for each cluster. 
From an obser vational perspective thes e results are broadly 
consistent with lDe Filippis et al.l (|201ll ) who compared what 
effect local and global background corrections have on the 
counts in clusters at low redshift. 

The radial distances for the mock sample are calculated 
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with respect to the central galaxy of the halo, which is then 
excluded from the fitting process. Once the best-fitting val- 
ues of K* are derived, we are able to evaluate JV500 (out to 
K — K* + 2) for each mock cluster and fit a power law of 
the form: 



1 00.0 F 



log N500 = a + P log M500 ■ 



(11) 

As described in Section \4. 2 1 the choice for using such a sim- 
ple relation for fitting the data points in our plot, is justified 
by the fact we analyse only satellite galaxies within haloes 
of mass (M200 > 10 14 M©) significantly higher than the cut- 
off halo mass (>, 10 12 Mp) introduced in the m odels (e.g., 
iKravtsov et~al]l2004l : IZheng. Coil fc Zehavi|[20(57h . 

In this way we are able to investigate the slope of 
the -/V500 — M500 relation and compare the obtained value 
from the simulations, with those obtained in the literature. 
This process is carried out for each of the aforementioned 
simulation snapshots, allowing us to study the evolution 
of the slope f3 of the -/V500 — M500 relation with cosmic 
time. The best-fit values of /3 are reported in Table [3] 

In order to investigate the evolution of the mean number 
of satellite galaxies (JV500} per halo mass with cosmic time, 
we restrict ourselves to the average values ({N500}', (A/500}) 
for the 100 most massive clusters in each snapshot and carry 
out the same process utilized for the observed sample, now 
based on the stacked LF of the mock clusters. The values of 
{N500) obtained are reported in Table [3] and plotted against 
M 500 in Fig. 



5.3 Concentration Parameter 

For each set of mock clusters, the value of c g is also calcu- 
lated using the same maximum-likelihood fitting procedure 
used for the observed clusters. The model used for the fit- 
ting procedure is similar to that described in Section [4.2.21 
but with Ebacfc = (see Eq. [9}, since the mock clusters 
do not suffer from contamination by galaxies not belong- 
ing to them due to the unambiguous way they are assigned 
to the mock haloes. The fitting procedure is applied to all 
the clusters individually (see Fig. [7|) and also applied to the 
stacked projected number density profiles, which are first 
normalized to the virial area within -R200 (see Table [3]). Al- 
though at a given halo mass, the concentration parameters 
are commonly assumed to be lognormally distributed (e.j 
|jinell200d ; iNeto et al. I l2007l ; IComerford fc Nataraianll2007 
we prefer in this work to use the average values of c g from 
the stacked profiles to compare with observations because 
the stacking procedure has the advantage of erasing individ- 
ual d eviations, usuall y related to the presence of substruc- 
ture (|Gao et al.ll2008T ). These deviations may contribute to 
the already considerable intrinsic scatter (|Neto et al.ll2007l 
treated this issue in details) in the parameters fitted to in- 
dividual profiles, possibly masking underlying trends in the 
data. It is worth pointing out that for the stacked profiles 
only the radial range 0.05 < R/R200 i$ 0.65 is used in the 
fitting process in order to avoid the innermost radial region 
where the NFW model fails to correctly reproduce the den- 
sity profile, as shown in several studies in t he literature (e.g ., 
iNavarro et "aT1l2004l ; iPrada et al1l2006l and iGao et aljfeoosf) . 
The best fit values of the average c g are reported in Table 
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Figure 7. (1 + z)c g against M2oo- Measures of c g are derived 
individually for the 100 most massive mock clusters in each of 
the considered snapshots. They are plotted as: blue squares (z = 
1.27), black squares (z = 0.99), grey squares (z = 0.83), red 
squares (2 = 0.62), magenta squares {z = 0.41), cyan squares 
(z = 0.21), green squares (z = 0). Average values are plotted as 
well: black full dots (mock clusters with z > 0.8), blue full dots 
(mock clusters with z < 0.8), grey full dot (mock clusters with 
2 = 0), green full dot (observed clusters). 



6 RESULTS 

6.1 Evolution of K* 

The K* values for our 15 clusters in the 3 redshift bins 
0.85, 0.97, 1.23 are shown in Table[2]and compared to other 
results in Figure U Our values of K* are consistent with 
other observations at similar wavelengths of other high- 



redsh i ft clusters (e.g., 
2004 IStrazzullo et al 



De Propris et al . 1999: El lis fc Jonesl 
120061 : lAndreonl l200d : iMuzzin et all 
2007al ). From a direct comparison with the K* values pre- 
di cted from ou r analy sis of the mock clusters in the SAMs 
ofH owcr et all (|2006l ). we find good agreement: indicating, 
for the first time, that the SAMs do a reasonable job of 
predicting K* evolution in rich clusters at redshifts z ^ 1. 
However, our data cannot distinguish between the mass- 
assembly SAMs and pure passively evolving SSP models 
with formation redshift Zf ~ 5, implying that red, old 
massive galaxies are already ubiquitous in rich clusters at 
z ~ 1.5. 



6.2 HOD 

6.2.1 Observations 

The observed (iVsoo) values measured from our clusters 
(0.8 < z < 1.5) are shown in Fig. \Sj where the lo cal data 
from 93 clusters from iLin. Mohr fc Stanford! (|2004 ) are also 
plotted along with their best-fit A^soo — M500 line. At all 
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Figure 8. Ratio of measured iVsoo over the predic i ton o f 
the low-z best-fitting relation of iLin. Mohr fc Stanford! j2004) , 
against z: observed clusters (black full dots), mock clus- 
ters with z > (red full dots), mock cluste rs with z = 
(green full d ot), ILin. Mohr fe Stanford! | |2004) data points 
(gr ey full dots), |Pe Pr opris e t al. I jl999h sample as reanlaysed 
by ILin, Mohr fc~sTam^rdri2004l) (empty squares) , iMuzzin et alj 
ll2007bh (crosses, note these values are extrapolated out to K* + 
3). Note our simulations based measures are average values mea- 
sured over the 100 most massive clusters in each snapshot. 



redshifts our points lie below the local fit. In Figure [H] we 
show the variation of (iVsoo) with redshift. Here the val- 
ues of JV500 have been normalised to the low-z best-fit re- 
lation shown in Figur e [5] using the relation showe d in Ta- 
ble 1 (second row) in lLin. Mohr fc Stanford] |2004| (see also 
Fig. 9 in their paper). The values measured for our clus- 
ters have 0.35 < A'soo/A'fit, local < 0.8. This result is in stark 
contrast to t he results of the re-analy sis by Lin and collabo- 
rators of the |Pe Propris et all (|l999l ) intermediate- 2 cluster 
sample, also shown in Figure [8] which indicate normalised 
■^500 values typically between 1 — 4 for clusters in the range 
0. 2 < z < 0.8. However, our results are consistent with those 

of m 

uzzin et al.l (|2007bl ) . who carried out a similar study 
to ours on 15 CNOC1 clusters at 0.19 < z < 0.55, and 
of lAndreon et a l. (2008), who found evidence of a possible 
break of the cluster scaling relations at 2 ~ 1. 

To quantify the dependence of -/V500 on z and cluster 
m ass in our data, w e fit the -/V500 data to the relation used 
bv lLin et afl (|2006h ■ and given by 



N(M, 



N (i + zy(M/M y 



(12) 

as in lLin et al.ll2006l ) are 



where N Q and M (set at lO 14 ^ M 
the normalization factors of the relation. In the fit we use 
our binned results for the 15 high redshift clusters and the 
93 low redshift data points of ILin. Mohr fc Stanford! 12004 
The best-fitting values are No = 53 ± 1, 7 = — 0.611q'2q and 
s = 0.86 ± 0.05. Figure [9] shows the contours corresponding 



Figure 9. Joint 68 and 95 per cent confidence regions for 7 and 
s (see eqn. 1121 . Black full point shows the best-fitting value for 
7 and s (-0.61, 0.86), while the cross stands for the center of 
the contours obtained during the marginalisation. Empty square 
represents the best- fitting values of 7 and s (plotted together with 
their l a errors) by I Lin et al.1 II2006T ) when including also the low- 
2 cluster sample used here, in addition to their intermediate-^ 
cluster sample. 



to the 68 and 95 per cent confidence regions for 7 and s. 
The quoted errors on 7 and s are calculated by projecting 
the 68 per cent contour onto 7 and s axes respectively. 

Comparab le measurements at our redshifts are those of 
ILin et al-lBoOfj who analysed 27 clusters with < z < 0.9. 
Our s value agrees reasonably well with their estimate 
(s ~ 0.8 ± 0.5) although evidence for an evolutionary trend 
in their data is significantly weaker (7 = —0.03 ± 0.27). 
In fact, our value of 7 indicates a significant trend with 
2 at ~ 3d level; although with only 3 points at z ^ 0.85 
more data and better redshift sampling are required to 
estab lish accur a te tre nds. We are in fact also consistent 
with ILin et all (|2006l ) at ~ 2a l evel. Unfo r tunate ly we 
cannot test our results by using ILin et al.l (|2006r i data 
points to repeat our fit, since their values are calculated 
within 7?20oo- 

Turning to the study of the concentration parameter, 
the value of c g = 2.8j] 8 found for the observed sam- 
ple at z ~ 1, is consisten t within la with the one of 



2.9otg:£) 



iLin, Mohr fc Stanford! (|2004T > at z ~ 0.06 (c g 
However, the latter value is quite different from the one s 
measured at similar redshifts (e.g.. iRines fc Diaferid [20061 ). 
In Figure [TD] the current measure of c g is plotted, together 
with the values found in the literature, as a function of 

log(2 + 1) - 1 ir— , 

A comparison with the predictions of lGao et al.l (| 20081 ) 

for dark matter haloes (see Section f6.2.2p at the same red- 
shift as the median of the clusters studied here is also carried 
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Figure 10. Average galaxy concentration parameter against 
log(z + 1): this work (black full dots for m ock clusters and 
blue full dot for our observed sample). iRines fe Diaferiol 
j2006h (red asterisks), iKatgert. Biviano fc Mazurd j2004h 
(blue asterisk, slightl y shifted towards lower z to improve 
data poi nt visibility). iBiviano fc Girard 3 d2003ll (black empty 
square). IBivia no fc P oggiantil l|20ld l (green em pty squares), 
iLin. Mohr fc Stanford! \2004 \ (black asterisk), IMuzzin et"afl 
l|2007ah (bla ck crosses, sl ig htly shifted towards higher z 
for clarity). [ Carlberg et al.l il997T ) (black empty rhombus), 
van der Marel et alj j200cf ) (black empty triangle). Note for 
mock clusters, the 100 most massive clusters in each snapshot 
are used for this measure. Trend lines for the simulation-based 
average values of c g measu red in this work ( black dashed line) 
and for the predictions of iGao et al ] l|2008h derived for dark 
matter haloes with the same halo masses as our mock clusters 
(red dashed line) arc also plotted. 



out, showing that our value of c g 
this prediction (cd m = 3.55) (see Fig. LH 



2.8_ () g is consistent with 



6.2.2 Observations & Simulations: Most Massive Clusters 

Turning to the results obtained for the simulations, the 
values of -/V500 measured from mock clusters (0 < 2 < 1.3) 
confirm the picture shown by the observed sample. In 
fact, as for the observed clusters, all the simulation-based 
values of -/V500 are lower than those predicted by th e 
low- 2 best-fitting relation of ILin. Mohr fe Stanford! (|2004 ). 
unlike those measured by Lin et al. for the sample of 
|Pe Propris et all (| 19991 ) (Figs. [5] & Table [3]). 

When the evolution of /3 with z, as determined from 
mock cluster samples (Fig. Ill|l . is studied, the values 
obtained seem to reveal no significant evolution with 
2 (as also indicated by a Spearmann's rank correlation 
test, which gives a correlation coefficient of r = — 0.2 
with a significance of its deviation from zero p = 0.6). 



Figure 11. Halo HOD power-law index, as measured from mock 
clusters in each snapshot, against (z+1). Note all the clusters 
with M200 > 10 14 Mg in each snapshot are used to carry out 
this measure. 



The value measured at z = j/3 = 0.820 ± . 008) i s 
consistent with the ones o f ILin. Mohr fc Stanford! (|2004 ). 
IMuzzin et ail (|2007bl ) and IPopesso et al.l (|2007t ) Despite 

many s tudies, e.g. 



this^ as noted in the introduction 



iMagliocchetti fc Porcianl (l2003t); IZehavil (|2004h ; IRozo et all 
|2007l ) and iPoggianti et al.l (|201of T m easure discordant /3 
value s ranging from j3 = 0.55 ± 0.043 |Marinoni fc Hudson! 
l2002f ) to /3 = 1.1± 0.09 dKochanek et al.ll2003l ). Part icularly 
high are the values found by Abbas et al.l (|2010l ). who, 
depending on the magnitude depth used and the redshift 
(0.1 < z < 1.3), found 0.99 ± 0.10 < /3 < 1.66 ± 0.18. 
However, the majority of the studies in the literature agree 
with j3 being inconsistent with unity (this might not be 
the case for some other simulation-based studies, but a 
meaningful comparison with these studies is particularly 
hard) . 

In order to compare the results obtained for cluster 
galaxies with dark matter, the predic tions of the Cd m (z, M) 
relation obtained bv lGao et all (|2008f l within the MS can be 
used as a reference. Gao and collaborators presented in their 
study best-fitting Cd m (A^20o) relations at z = 0,0.5,1,2 
and 3. To carry out this comparison, the values of Cd m 
predicted by their model are calculated using the values 
of M200 at 2 = 0, 0.62 and 0.99 for our MS mock clusters 
(Table [3] Fig. 1101) and at 2 ~ 1 for our observed cluster 
sample (cd m = 3.55). This makes it possible to calculate 
Alogc = log(cg/cdm) (|Zheng fc Weinberg! 120071 ) at each 
of the chosen z (Table [3]) and to compare galaxy spatial 
distribution with the one of dark matter (a positive value of 
A log c would mean galaxies are more clustered than dark 
matter and vice versa for negative values). This test shows 
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that Alogc decreases with redshift, reaching negative 
(Alogc = —0.1, observed sample) or null (Alogc = 0, 
mock clusters) values at z ~ 1. At this z, the value of 
A log c is < 5 per cent of its value at z = 0. 



6.2.3 Simulations: The Effect of Mass Selection 

The results for ./V500 and c g shown in the previous sections 
are acquired using mock samples which have an increasing 
mean cluster mass with cosmic time (see Table [3]). As 
a result it is possible that this mass evolution affects 
the evolution of the HOD with cosmic time. In order to 
investigate this we repeat our analyses on two further mock 
samples of 100 clusters per snapshot: (i) selected within a 
constant mass range over cosmic time, set by the high-z 
sample; (ii) selected to have the same average mass at z = 1 
and z = 0. We discuss the results from these samples below. 

(i) We repeat our analysis on 100 mock clusters per 
snapshot randomly selected within the mass range of our 
observed cluster sample (0.4 < Msoo < 6.2 x 10 14 M Q ). 
The values of the average masses of these mock cluster 
samples are M 500 = 1.58, 1.45, 1.34, 1.21, 1.12, 1.06 and 
1.04 x 10 14 M Q at z = 0, 0.21, 0.41, 0.62, 0.83, 0.99 and 
1.27, respectively. When using these new mock clusters, 
we find very similar results to those described in Section 
16.2.21 for the 100 most massive mock clusters at each 
snapshot. In particular we find N^oo/Nat, local is always 
< 1 and for the 7 previously mentioned snapshots we 
find: iVsoo/AfiMocai = 0.39 ± 0.02, 0.42 ± 0.02, 0.45 ± 0.02, 
0.49 ± 0.03, 0.35 ± 0.03, 0.37 ± 0.03 and 0.53 ± 0.03 
respectively. 

We also recalculate the average concentration parameter 
using the 2 = and z = 1.27 snapshots and again reproduce 



similar results within the errors: 
and c g = 8.7±^3 at z = 0. 



4.7 



+0.97 



at z 



1.27 



(ii) As an additional test, we select mock clusters in 
the two redshift snapshots corresponding to z — and 
2 = 1.0. At both these redshifts we generate a sample of 
100 mock clusters, selected at random from a Gaussian 
centered at M500 = 2.1 x 1O 14 M0, with a width given 
by the dispersion of our observed high-redshift clusters. 
Recomputing values in the same way as described in 
Section RT2~21 we find (A 50 o) z=1 = 28 ± 1 at z = 1 and 
(A500) n = 36 ± 1 at z = 0. Norma lising to the low-2 rela- 
tion of iLin. Mohr fc Stanford! (|2004T l as before (see Fig. [S]), 
we find again values less than 1: iVsoo/^Vfit.iocai = 0.62 ±0.02 
and A 5 oo/A flt ,iocai = 0.5 ± 0.02 at z = and 2 = 1 
respectively. 

For the concentration parameter we get c g = 7.5^Q'g5 and 
c g = 3.7^0 44 at z = and z = 1, respectively. These results 
lie no more than 1 — 2a from the values at the corresponding 
redshift snapshots shown in Fig. 1101 while the larger errors 
reflect the increase in Poisson noise due to the smaller 
number of galaxies in less rich clusters, which are more 
numerous in samples characterised by lower average masses. 

The results obtained here for -/V500 and c g are both 
consistent with the simulation-based measurements in Fig. [8] 
and Fig. [10] respectively, which were based on mock samples 
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Figure 12. Mass distributions of lLin. Mohr fc Stanford! 12004 ) 

cluster sample (full line) and its mass matched counterpart in the 
mocks (dashed line, offset of 0.1 dex toward higher masses for 
clarity). See text for a description of the latter sample. 



whose average mass changes with epoch. They demonstrate 
the robustness of our claim that high-2 mock clusters agree 
with observations in hosting smaller numbers of galaxies 
compared to previous results and the claim of a significant 
evolutionary trend for c g . This conclusion is supported 
by other evidence in the literature which shows that the 
dependence of the concentration paramet er on mass i s wea k 
for haloes more massive than 1O 14 M0, e.g. lGao et al.l (|2008T ). 

Finally, we note that the HOD power-law index ft evo- 
lutionary analysis shown in Fig. [TTJ uses all clusters above 
M200 = 10 14 Mq, in order to maximise the mass range and 
increase the signal given the large intrinsic scatter of the 
N— M relation, with the consequence that the average mass 
in any of the 7 redshift snapshots varies by less than 30 per 
cent from the overall average across all snapshots. 



6.2.4 Simulations: A Direct Comparison with Lin et al. at 
z = 

The mock sample at z = formed by the most mas- 
sive haloes has a mass range (5.2 < M500 < 26.0 x 
10 14 Mq) which onl y part ially overlaps with that of 
ILin. Mohr fc Stanford! i|2004h (0.2 < M 500 < 12.6 x 
10 14 M Q ), resulting in a lack of haloes in the mocks less 
massive than ~ 5 x 10 14 M©. To investigate whether this 
disparity in the mass range is related to the large differ- 
en ce we find between our 2 = mock value of c g and that 
of ILin. Mohr fc Stanford! (|2004 ). we generate 1000 clusters 
selected randomly from Lin's accumulated mass distribu- 
tion and associate each of them with the 2 = mock clus- 
ter closest in mas s. This produces a mock clu ster sample 
closely matched to lLin. Mohr fc Stanford! (|2004h ; see Figure 
1121 When using the new mock sample described we obtain 
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c g = 8.01q3i, which is again consistent with the previous 
results for the low-z mocks and significantly different fr om 
the value (c g = 2.9±o,2a) of iLin. Mohr fc Stanford! <j2004l ). 



7 DISCUSSION 

From our study of the evolution of K* , shown in Fig. [4] it 
is evident that the SAMs reproduce the K* values of our 
high redshift clusters reasonably well. The data points of 
ICirasuolo et al.l (l2007h and the paramet erized fit of the evo- 
lution of K* from lCirasuolo et al.l (|20ld ) for the field LF are 
also plotted in this Figure. It can be seen there is a small 
difference of, on average, AK* ~ 0.5 at z > 0.8 between 
the cluste r environment (brig hter K*) and the field. Fur- 
thermore, iBalogh et al.l (|200lh found in their study of the 
J-band LF a difference of similar magnitude between clus- 
ter and field of AJ* ~ 0.5. The results of Cirasuolo et al., 
are from a sample of 50,000 galaxies with photometric red- 
shifts 0.25 < z sC 3.25 from the UKIDSS UDS in which they 
compar e the jf-band LF w ith the predictions of SAMs (in- 
cluding iBower et al] |2006h finding broad consistency - even 
though some models still tend to under predict the num- 
ber of field galaxies at high z. However, we point out that 
the I Bower et al. I (|2006h SAM does a better job of predict- 
ing the number of high-z field galaxies than other SAMs. 
Although our study is restricted to the analysis of K* , it is 
clear that our results indicate that the current SAMs can 
predict the luminosities of L* galaxies both in the field and 
rich clusters, where substantially more merging has taken 
place , reasonably well. On the other hand some author s 
(e.g., ICimatti. Daddi fc Renzinll200d : IPozzetti et"afll20ld ) 
have pointed out that the massive end of the galaxy mass 
function does not appear to evolve significantly since z ~ 1 
and questioned whether downsizing is compatible with the 
predictions of SAMs. As mentioned in the introduction, the 
situation is potentially even worse for BCGs at the centre 
of rich clusters. In fact, Collins et al (2009) and Stott et al 
(2010) showed that the observed stellar mass of BCGs in 
high-mass X-ray clusters does not evolve significantly since 
2 ~ 1.5, in contrast with the SAMs which predict that BCGs 
at this redshift should only have 20—30 per cent of the stellar 
mass of their local counterparts. 

These results, taken together with the small difference 
between observe d field and cluster K* values, now confirmed 
here, (see also I Cole et all l200ll ; |Pe Propris fc Christleinl 
l200d ) and together with other evidence found in the litera- 
ture that the luminous (L > L*) part of the LF weakly de- 
pends on the environment ( e.g., ICimatti. Daddi fc Renzinl 
l200d ; iDe Filippis et alj201ll ) , indicates that the hierarchical 
galaxy formation scenario works well for galaxies up to a 
mass M ~ I0 11 Mq (about the mass of an L* galaxy). 

Our halo HOD power-law indexes (/3 ~ 0.8) for 
mock clusters are in agreement with s evera l studies in the 
literatu re, e.g. ILin. Mohr fc Stanford! |2004h , iMuzzin et al.l 
(|2007bh and IPopesso et al.l (|2007fi . However, the pic- 
ture is complic a ted, since several other studies such as 
iKochanek et all (|2003h and IPoggianti et all (|2010D . found 
values of /3 greater than or consistent with 1. As emphasized 
by Popesso and collaborators and Lin et al., a value P < 1 
is expected from the theoretical point of view. This is 



because, despite the fac t that the hierarchic al structure 
form ation models (e.g. IDe Lucia et alj |2004| ; iGao et all 
l2004T l predict a universal mass distribution of sub-haloes 
(independent of the parent halo's mass and with the con- 
sequence that the number of sub-haloes is proportional to 
halo mass), the introduction of baryons produ ces a decreas 



ing n umber of galaxies per unit mass (e .g., Berlind et al 



2003) and an increasing M/L ratio (e.g., iKauffmann et al 



1 9991 ) for higher halo masses. Several processes may be 



respo nsible for this behavior, such as an increased merger 
rate l|White. Hernquist fc Springell 200lf) or an i ncrea sing 
galaxy destruction rate (|Lin. Mohr fc Stanford! I2004T ) or 
decre asing star formatio n and gas cooling efficiencies 



(e.g., Berlind et al 



20031 ). However, as pointed out by 



IPopesso et alj (2007), each of these processes should 



leave their mark in t he p roperties of galaxy clusters. 
ILin. Mohr fc Stanford! (J2004) split their sample into a 
high and low-mass sub-samples and found both an excess 
and a lack of galaxies in the faint and the bright parts of 
the LF respectively for the low-mass cluster sub-sample, 
although they concluded that none of the studied processes 
(tidal stripping, ram pressure stripping, galaxy harassment, 
variations in star formation efficiency, stellar aging, dif- 
ferences in galaxy radial distribution, dynamical friction 
and galaxy mergers) can by itself satisfactorily explain 
these differences. On the other hand, IPopesso" et all (|2007f) 
found a universal LF with no dependence on the cluster 
mass and wave-band. For this reason, they proposed the 
only way to have P < 1 leaving the properties of galaxy 
clusters unchanged is if the sub-haloes mass distribution is 
not universal. It is beyond the scope of the present study 
to address which process may cause /3 < 1 or whether the 
sub-haloes mass distribution is universal and we simply 
point out that our analysis on mock clusters shows no 
significant evolution of /3 with z (/3 ~ 0.8 out to z ~ 1.3) 
and that whatever this process might be, it is already in 
place by z ~ I . 

Our observational-based results indicate that low-z 
clusters are richer than their high-z counterparts when com- 
paring clusters of the same mass. Using the best fit relation 
to Eq. II2l and the masses in Table O iVsoo increase by a fac- 
tor of ~ 1.6 from z ~ 1 to z ~ 0.6 and by a factor of ~ 2.6 
from z ~ 0.6 to z = 0. This growth rate of the number 
of gal axies is very similar to the one found by I Abbas et alj 
|20ld ) over a similar redshift interval. 

Using mock clusters, IPoggianti et~aH |20ld ) also found 
little evolution in the /3 value for their samples, although 
pointed out that values can depend on whether masses of 
mock clusters are derived from their intrinsic halo mass 
(P = 1. 00 ± 0.04) or using their velocity dispersion mass 
proxy (/3 = 0.77 ± 0.03). It is not possible to match our 
X-ray measurements to the simulations in the same way 
at the present time but clearly further work is necessary 
to compare different mass estimators. However, we point 
out that we obtain an average value of P < 1 for our mock 
clusters despite using the halo mass values provided by the 
simulation. 

Focusing the attention on the galaxy concentra- 
tion parameter (Table [3] and Figs. [6] & I10|) . the present 
analysis leads to a picture which is difficult to fully un- 
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large scatter is probably due to the heterogeneity of the 
analysis methods, which are often based on different radial 
distances and magnitude limits. Another possible cause of 
this large scatter could be the cluster-to-cluster variation, 
which can influence the headline value of c g when simply 
averaging over individual measurements (using individual 
surface density profiles), as opposed to the preferred method 
of fitting the stacked density profile. However, we note 
that the most accurate c g estimate o n intermediat e z scale 
(c g = 4.13 ± 0.57), carried out by iMuzzin et all (|2007aT ) 
and based on redshift measurements, agrees well with our 
observations at z = 1. The values of c g measured for our 
mock clusters are seen to increase with decreasing redshift, 
giving rise to a significant correlation with a Spearman 
rank correlation coefficient of —0.96. In particular, those 
obtained for high-z mock clusters (c g = 3.851q'4 5 , 3.70lo'| 4 
and 3.77±g;|f at z =0.83, 0.99 and 1.27, respectively) are 
found consistent within la with the value of c g = 2.81J g 
found for our observed sample at z ~ 1. The large errors 
characterising the observed values found in the literature, 
do not allow us to address this issue in a more quanti- 
tative way. Having said that, when the value obtained 
for mock clusters at z — (c g = 6. 301°'^) is compared 
with the observed one of llin. Mohr fc Stanford! l|2004 ) 
(c g = 2.90^0 22) a t 2 ~ 0.06, the discrepancy appears large 
and significant, although the models agree with other low-2 
observation-based estimates of c g . If real, this would be 
an indication that SAMs have problems in predicting the 
evolution of cluster galaxies observed from z ~ 1 to z — 0. 



A potentially very important difficulty which arises 
in the simulations is the problem of the orphan galaxies 
(Shaun Cole, private communication). These are galaxies 
which have lost their dark matter haloes and whose tracks 
in the simulations are traced by following the particle that 
was the one most bou nd in its halo prior to disruption. In 
the lBower et all (|2006h SAM, the merging of satellites onto 
the central galaxy is treated using a dynamical friction 
calculation that, instead of using the true orbit of the 
galaxy, uses random orbital parameters that are only 
statistically consistent with N-body simulations. This could 
lead to biasing the concentration parameter towards higher 
values. In fact, in reality, the satellites most likely to merge 
are the ones on orbits closest to the centre, whereas random 
selection in the SAM of Bower et al. removes satellites 
randomly at all radii. Hence, the remaining satellites may 
be too concentrated. While there is no reason to expect that 
the current treatment of satellite galaxy mergers should 
give rise to the entire trend seen in Fig. 1101 the possibility 
cannot be ruled out. This issue will only be decided by 
carrying out a more accurate treatment of the merging 
and orbits of satellite galaxies, which has already been 
developed, but not yet applied to the MS. However, it is 
worth pointing out that the problems of cal culating the 
radial distribution of satellite galaxies in the iBower et all 
|2006l ) model are not expected to lead to significant errors 



in the numbers of satellite galaxies contributing to the HOD. 

The growth of mock c g with decreasing z seen in 
the present analysis is not unexpected. What is sur- 
prising, instead, is the rate of this growth. Investigating 
dark matt e r ha lo es within simu lat ions, studi e s like 
iNetoetalJ (I2007T); iGao et all J2008I ); iDuffv et all (|2008h 
and iMunoz-Cuartas et all (|201lT ) addressed the Cd m (z,M) 
relation in details finding Cd m being anti-correlated with 
mass and z. In particular, the anti-correlation between Cdm 
and M weakens at higher z and evidence of dark matter 
haloes h aving a con s tant C dm ~ 3.5 — 4 at z > 1 have been 
found bv lGao et all (|2008h . 

Comparing the mock-cluster values of c g measured here 
with the ones o f the dark matte r haloes predicted by theo- 
retical studies (|Gao et alj|2008h at different z (see Section 
O, the results seem to show that galaxies have a similar 
concentration to dark matter at z ~ 1 but subsequently 
they become more concentrated with decreasing redshift. 
This may be due to the effect of non-gravitational or 
gravitational processes, such as: gas cooling, AGN feedback, 
dynamical friction, merging and tidal stripping, which have 
significantly modified the galaxy distribution within the 
dark matter halo over the last ~ 9 Gyr. This growth is 
domi nated by the cont inuously refined physics of the SAMs 
(e.g., IDuffv et alJboid l and given this sensitivity, a glance 
at the data in Fig. [10] suggests the need for a thorough and 
homogeneous study of c g in observed galaxy clusters and 
groups over a wide range in redshift. 



8 CONCLUSIONS 

We study the evolution of cluster galaxies with z using K- 
band photometry of a cluster sample made of 15 of the 
highest-z (0.8 < z < 1.5) X-ray clusters observed so far. 
In particular, we investigate the K* Hubble diagram out to 
z ~ 1.3 in comparison with SSP (no evolution in absolute 
_R"-band magnitude and passive evolution) and SAM evo- 
lutionary models. This allows us to explore the process of 
galaxy formation. We also study the HOD of this sample, 
by investigating (N), as a function of M and z, and c g with 
z in comparison with Cdm- All our analysis is carried out in 
strict comparison with simulations, through the use of mock 
cluster samples taken from MS at < z < 1.3 and populate d 
with galaxies by means of the SAM byH OWCT et all |2006h . 
In addition to being an important test for SAMs, this allows 
us to investigate in details the evolution of f3, (N) and c g 
with z and to study how the relationship between galaxies 
and dark matter spatial distributions within haloes changes 
with cosmic time. This kind of study is fundamental for un- 
derstanding the role non-gravitational processes, as part of 
the physics of galaxy formation, take in influencing the way 
galaxies populate dark matter haloes. Our main conclusions 



(i) The values of K* obtained show the existence of old, 
evolved and massive galaxies at z > 1.5. Despite the evolu- 
tion of K* is well reproduced by SSP passive evolutionary 
models (zf > 3), we can not disentangle b etwee n these 
models and the predictions of IBower et all l|2006h SAM, 
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which is seen, for the first time, to reproduce well the 
evolution of K* also in the cluster environment. 

(ii) When comparing the values of K* obtained for ob- 
served clusters with the ones of the field, there are no 
major differences. This fact, coupled with the marginal 
evolution of the massive end of the galaxy mass function 
seen observationally (e.g.. Icimatti. Daddi fe Renzinil 120061 : 
IPozzetti et alj|20ich . leads to questioning whether the role 
of the environment is negligible for the formation of galaxies 
more massive than ~ 10 11 Mq (about the mass of a L* 
galaxy). 

(iii) By investigating A500 in real and mock clusters, 
this study shows that high-z clus t ers a re poorer than 
those at low z. Unlike iLin et all (|2006T l. we find sig- 
nificant trends of A500 with both z and cluster mass: 
N(M,z) = (53±l)(l + 2)"°- 61 -^o(M/10 14 - 3 ) a86±0 06 . 

(iv) Using mock clusters, the slope /3 of the -/V500 — A/500 
relation is found to be significantly lower than one at all z 
(out to z ~ 1.3), showing no significant signs of evolution. 
This means more massive clusters are characterized by a 
lower galaxy number per unit mass compared to lower mass 
systems already at high z. Because of this, the local value 
of P < 1 can not be explained as being due to local clusters 
found to be richer than those at high-z. 

(v) Although our results here and those from the literature 
seem to indicate a decreasing trend of c g with z, overall 
the data and SAM-based predictions are very uncertain 
and the current situation emphasizes the compelling need 
of a systematic and homogeneous study of the galaxy con- 
centration parameter in clusters over low and intermediate 
redshifts. 

(vi) When comparing our mock-cluster values of c g (3.7-6.3) 
with those of Cd m (3.6-3.8), starting from similar concen- 
trations at z ~ 1, galaxies seem to become more concen- 
trated than dark matter as z approaches 0. This may be 
due to gravitational and/or non-gravitational processes sig- 
nificantly modifying the distribution of galaxies within dark 
matter haloes over the last ~ 9 Gyr. However, the problem 
of the orphan galaxies in the SAMs prevents us concluding 
that this is a real physical trend. 
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